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ABSTRACT 

Proplyds are photodissociation region (PDR)-like cometary cocoons around young stars which are thought 
to originate through photo-evaporation of the central protoplanetary disk by external UV radiation from the 
nearby OB stars. This letter presents spatially resolved mid-infrared imaging and spectroscopy of the proplyd 
HSTIO obtained with the VLT/VISIR instrument. These observations allow us to detect Polycyclic Aromatic 
Hydrocarbons (PAH) emission in the proplyd photodissociation region and to study the general properties of 
PAHs in proplyds for the first time. We find that PAHs in HSTIO are mostly neutral and at least 50 times 
less abundant than typical values found for the diffuse ISM or the nearby Orion Bar With such a low PAH 
abundance, photoelectric heating is significantly reduced. If this low abundance pertains also to the original 
disk material, gas heating rates could be too low to efficiently drive photoevaporation unless other processes 
can be identified. Alternatively, the model behind the formation of proplyds as evaporating disks may have to 
be revised. 

Subject headings: circumstellar matter — ISM: lines and bands — protoplanetary disks — stars; individual 
(HSTIO) — stars: winds, outflows 



1. INTRODUCTION 

Most lo w mass stars are b om in transient OB associa- 
tions (e.g. JLada & L ada'2003). and there is evidence that our 
own young Solar System evolved near massive stars (e.g., 
iHester et al.l 120041). Externally illum inated protoplanetary 
disks or proplyds ( lO'Dell et alj|1993l) are young stellar ob- 
jects (YSOs) surrounded by Solar System-sized protoplane- 
tary disks and found embedded within or near a H ii region. 
In these extreme environments the disks are exposed to in- 
tense UV radiation fields and stellar winds from the OB stars, 
dynamical encounters with sibling stars and supemovae, on 
timescales of planetary system formation and early evolution. 
Hence, the study of their properties may bring key constraints 
on the understanding of the general mechanism of planet for- 
mation and the origins of our Solar System. 

Proplyd morphology has been explained by models of evap- 
orating flows in externally illuminated disks or globules (e.g. . 
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FUV (6 eV< E < 13.6 eV) heated photodissociation region 
(PDR) encased within a EUV (E > 13.6 eV) ionized hydro- 
gen cocoon with a "head-tail" shape. FUV photons of nearby 
OB stars penetrate deep into the disk and, at relative high col- 
umn densities, dissociate hydrogen molecules, ionize carbon, 
and heat the gas to T~400-4000 K, forming the PDR. The 
resulting thermal pressure drives an expanding hydrodynam- 
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ical flow of neutral material escaping from the disk surface 
with velocities of 1-3 kms"'. The mass-loss rate generated 
by this photo-evaporation wind determines the lifetime of the 
gaseous disk and hence the timescale for the formation of gi- 
ant planets. At some distance from the disk, the EUV pho- 
tons ionize the neutral wind and form an ionization front (IF) 
with T~10"* K. The observed tails result from the diffuse UV 
radiation field (produced through the recombination of sur- 
rounding nebular gas) which drives an evaporation flow on 
the shadowed side of the disk. These models consider the 
FUV photoelectric effect on sm all dust grains and Poly cyclic 
Aromatic Hydrocarbons (PAHs. jjoblin & Tielensl201 ll) as the 
main gas-heating mechanism at the disk surface. Recent ob- 
servations (Okada et al., submitted) show that PAHs play 
a major role in the heating process, and theoretical studies 
jKamp & Dullem ond 2004) confirm this also for protoplane- 
tary disks. The proplyd models mentioned above are adapted 
from classical PDR models (see [Hol lenbach & Tielens 1997| 
for a review) which consider standard interstellar medium 
(ISM) PAH abundances. However, up to now, there is no ob- 
servational abundance determination for prop lyds. 

HSTIO (182-413; lO'Dell & WenI [19941) is a teardrop- 
shaped proplyd {I'.'l X 2'.'6) containing a prominent nearly 
edge-on disk (i~80°, PA~86°) visible as a dark silhouette in 
Ha, optical ionized species and the continuum, but glowing 
in [OI1 /16300 and in the 2.12 ^m ro-vibrational line of H2 
(iBallv et alJl2000l: iChen et al J [19981 Vicente et al., in prep.) 
This object is located at a projected distance of 56" from 0^ 
On C to the SE, the main ionizing star of the Trapezium Clus- 
ter (414+7 pc, Menten_eL.al, ,20Q7)- In this letter we present 
spatially resolved VLT/VISIR observations of PAH emission 
in HSTIO. We determine the PAH abundance in the neutral 
flow and discuss the implications for photo-evaporation mod- 
els. 

2. OBSERVATIONS AND DATA REDUCTION 

Mid-IR imaging (75 mas/pix) and low-resolution spec- 
troscopy (R~200-350, 127 mas/pix) of the proplyd HSTIO 
were obtaine d on the 13th December 2005 with VISIR 
jLagage et al] l2004h . under good ambient conditions (see- 



2 



Vicente, Berne, Tielens et al. 



Ha / H2 

/(J^ Trapezium ^ 

1" 1" 

HST/WFPC2 VLT/NACO: 




Figure 1. HSTIO images (from left to right): Ho- (0.656 pm) from HSTAVFPC2 (45.6 mas/pix; HST arcliive), H2 (2.12 /jm) from VLT/NACO (27 mas/pix; 
Vicente et al., in prep.), and PAHl (8.6 pm) and PAH2 (1 1.25 pm) from VLT/VISIR (75 mas/pix; this work). The VISIR images show the I'/O slit position for the 
8.6 and 1 1 .4 pm spectral settings respectively with the intensity scale given in mjy/arcsec^ . The arrow indicates the direction towards 6' Ori C, the main ionizing 
O-star of the Trapezium cluster North is up and east to the left. 



ing at 0.5 /^m between 0'.'6 and 1"). A spectro-photometric 
standard (HD 35536) of similar airmass was observed imme- 
diately before and after the science target to allow for tel- 
luric absorption correction and calibration. The data were re- 
duced with a VISIR customary pipeline based on IDL scripts 
(iPantin et al.ll2005|) and corrected for the background using a 
multi-resolution inpainting scheme (Pantin 2010). The im- 
ages were collected with filters PAHl (8.59+0.21 /jm), SIV 
(10.49+0.08 iim), and PAH2 (11.25+0.3 fim) for a total expo- 
sure time (on source) of 30min and using the parallel chop- 
ping/nodding jitter mode with a chopping throw of 7". The 
spatial resolution of the final images (19" x 19"), measured in 
the only point source in the field, was 0'.'35 for PAHl and 0'.'34 
for PAH2, the latter being diffraction limited and reflecting 
the improvement in the seeing HSTIO is spatially resolved 
in the PAHl and PAH2 images (a.'97x 2'.'5), but not detected 
in SIV after nebular subtraction. No central star inside the 
disk is visible in any of the VISIR images (Fig. 1). 

The A^-band spectroscopy consisted of three settings cen- 
tered at 8.5, 9.8 and 1 1.4 /im (Ih integration time) chosen for 
covering the lOjum silicate feature and PAH emission bands at 
8.6 and 11.25 fim. The 1" slit was placed along the proplyd 
head-tail, in the NS orientation, and the chop/nod was per- 
formed along the slit with a throw of 8". The ID-spectra were 
extracted by integrating the flux for each wavelength over the 
spatial extension of the proplyd (~17-19 pix) in the 2D back- 
ground subtracted spectra. They show spatially resolved PAH 
emission at 1 1 .25 fim (Fig. 2) but no detection above the noise 
level at 8.6 jum. The slit position, estimated using the RA and 
DEC values from the headers, is overlaid on top of the VISIR 
images in Fig. 1 and have an uncertainty of 0'.'2-0'.'3, corre- 
sponding to the accuracy of (small) relative offsets with the 
VLT. The RA offset of t he slit at 8.6Aim, relat ive to the pro- 
plyd nominal position in lO'Dell & Weill (11994 . is 0'.'6-a.'8 to 
the east which leaves the proplyd out of the slit. We believe 
this was due to errors occurring during the blind offsetting per- 
formed to obtain the HSTIO spectra. For the 1 1.4 fim setting 
the proplyd falls almost entirely inside the slit (see Fig. 1). 

The H2 2.12 fim image (27 mas/pix) was collected with 
the VLT adaptive optics instrument NACO under program 
ID 076.C-0874 (PI, S. Vicente). This data is described and 
analyzed in a forthcoming paper (Vicente et al., in prep.). The 
Ha 0.656/im image (45.6 mas/pix) from the Hubble Space 

^ The diffraction limit at the VLT (8.2m) varies from ^'26 to 0'/34 in 
the range 8.6 - 1 1.25,uiti. As the size of a UT miiTor is comparable to the 
turbulence outer scale. VISIR data are already diffraction limited for optical 
seeing below 0'.'6. 
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Figure 2. VLT/VISIR low-resolution spectrum (R~200) of the proplyd 
HSTIO showing the broad PAH band at 11.25 pm, and extracted from the 
2D spectruin on the top. The PAH2 filter range is indicated by the vertical 
dotted Unes. A flux of 16 mjy for the PAH plateau was determined as the 
mean flux in an interval of the same length as the PAH2 filter and is repre- 
sented by the horizontal dashed line. Over-plotted is the spectrum smoothed 
with a boxcar = 8. 

Telescope instrument WFPC2 was retrieved from the ESO 
archive (program GO 6603, J. Bally). The optical and near-IR 
images were rebinned to the same pixel size of the VISIR 
images and aligned with sub-pixel precision (0.01 pix) with 
IMALIGN/IRAF. The NACO image was used as reference 
because it contained several point sources in the field, con- 
trary to the VISIR and HST images for which only one to two 
stars could be found. Although there is a perfect overlap of 
the centroids of the stars in the RGB image, the estimated er- 
ror in the alignment is ± Ipix (or ± 75 mas) from comparison 
of the position of the silhouette disk seen in the Ha image to 
the H2 disk emission seen in the NACO image (Fig. 3). This 
error results from the different Point Spread Functions (PSFs) 
observed at the different wavelengths (0'.'06 for Ha, 0'.'08 for 
H2, and 0'.'35 for PAHs). 

3. OBSERVATIONAL RESULTS 

Fig. 1 shows the optical, near-IR and mid-IR images of 
HSTIO tracing the different key elements of its morphology. 
The optical Ha image traces the ionized gas at the ioniza- 
tion front with its peak surface brightness facing the Trapez- 
ium stars. This is a clear evidence that 0' Ori C is the main 
source of UV radiation driving the photo-ionization and evap- 
oration in HSTIO. The H2 2.12 fim emission traces the molec- 
ular gas at the disk surface. This was first discovered in a 
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Figure 3. Color composite image of the proplyd HSTIO applying a square-root stretch in intensities for better contrast, and showing the different regions traced 
by the different wavelengths: Ha tracing the ionized gas at the ionization front (blue), H2 the molecular gas at the surface of the nearly edge-on protoplanetary 
disk (green), and the PAHs emission the atomic gas within the PDR (red). From left to right we increase progressively the intensity of the PAH2 emission and 
find it to be coincident with the more extincted or less emitting regions in the Ha optical image. The silhouette disk is indicated by the green ellipse. 

The PAH emission is locaUzed within the ionized envelope 
and at the disk surface from where, according to PDR mod- 
els, the FUV-heated material evaporates generating the neu- 
tral wind and creating the proplyd PDR. Additionally, even 
thought the PAHl and PAH2 emission show similar head- 
to-tail extent, their distribution is different. PAHl emission 
at 8.6 fim is fainter (~ 2x in surface brightness) and more 
homogeneously distributed within the cocoon, whereas the 
PAH2 emission at 1 1 .25 jim is sharper and brighter on the op- 
posite side to the direction of 0' Ori C, and coincident with 
the areas showing less emission or more heavily extincted 
in the optical Ha image (Fig. 3). Considering that most of 
the optical extinction is caused by ~0.1yum dust grains, and 
these are expec ted to be depleted in th e neutral flow (with 
Av = 0.1 - 0.2. lHennev & 0'Delll[T999l) due to grain growth 
and settling in the protoplanetary disk, the 11.25 fim emis- 
sion is tracing the regions of higher density in atomic gas and 
possibly very small particles of dust. Additionally, while the 
11. 25 fim feature is dominant for neutral PAH molecules, the 
8.6 fim feature is stronger f or ionized PAHs raostly present in 
high UV irradiated regions (lJobUn&Tielensll2011l and refer- 
ences there in). Hence, the difference in spatial distribution 
observed in the two filters may reflect a charge effect associ- 
ated to the proplyd location in the foreground of the nebula, 
that is, in between 0^ Ori C and the observer Positively ion- 
ized PAHs are expected to be more abundant on the irradi- 
ated side of the proplyd opposite to us, whereas the bulk of 
the PAHs reservoir, as seen in the proplyd "shadowed" side, 
seem to be in the form of neutral molecules. The h.e/Iii.is 
band ratio, at the position where we extracted brightnesses in 
HSTIO, is of th e order of 0.5, in ac cordance to the astronom- 
ical template of iPilleri et al.l (120121) for which a value around 
0.4 is found for neutral PAHs and 1.45 for PAH+. For high 
radiation fields, PAHs can be neutral if the density of the gas 
in the flow is high (~ 10*" cm'-^. jTielens 2005), allowing for 
efficient recombination of PAH cations with slow electrons. 
As we will see, this is most likely the case (Sect. l4~2l i. Finally, 
we note that there is no PAH emission beyond the ionization 
front, suggesting that they are largely destroyed beyond this 
point as seen in the Orion Bar (Giard et al... 19941) . 



HST/NICMOS image ( iChen et al.|[l998l) and more recently 
confirmed with adaptive optics ground-based imaging (Vi- 
cente et al., in prep.). The VISIR mid-IR images reported here 
show emission associated to the proplyd in both filters PAHl 
(8.6 fim) and PAH2 (11.25 fim). The latter detection is also 
confirmed with spectroscopic observations (Fig. 2) and can 
be attributed to (solo) C-H out-of-plane bending mode in PAH 
molecules with long straight edges. No continuum is observed 
below llyum, but beyond 11.3 fim we see in Fig. 2 the 11-14 
fim PAH plateau resulting from the blend of C-H out-of-plane 
bending modes (e.g. arising from PAH clusters). On the basis 
of the similarity of the morphologies in the PAHl and PAH2 
images, and since the PAHl filter is also situated on top of 
a C-H vibration (in-plane bending mode), we conclude the 
PAHl emission is also due to PAHs. The total density flux in- 
tegrated over the full proplyd extension (rectangular aperture 
of 25 x45 pix) in the PAHl image is 26.6±0.5 mJy (3cr) when 
considering only the photon noise. Given that the uncertainty 
in the conversion factor used for calibrating the images is typi- 
cally 10%, we obtain 26.6+2.6 mJy. The integrated flux of the 
1 1.25 fim feature in the spectrum, divided by the passband of 
the PAH2 imaging filter (10.95-1 1.55 fim), is 29.4+12.9 mJy 
(Icr) which is consistent, within the error, with the extracted 
photometry in the PAH2 image, 3 8. 8 ±3. 8 mJy. The smaller 
value obtained from the spectrum with respect to imaging is 
likely due to slit losses which were worsen by pointing in- 
accuracies from the blind offsetting. The 11.25 fim emission 
observed in the 2D-spectrum in Fig. 2 is extended over ~17- 
19 pixels or 2'.'2-2'.'4 (127 mas/pix), similar to the head-to-tail 
size of HSTIO in the PAH images, and is about 50 mJy at 
the peak position. The maximum surface brightness in the 
PAHl and PAH2 filters (/g.e and /i 1.25) were derived by taking 
the average of the four brightest pixels at the proplyd head in 
each VISIR image. These values are given in Table 1. 

4. PROPERTIES OF PAHS IN HSTIO 

4.1. Spatial distribution and ionization 

In PDRs, PAHs emit mostly in the neutral gas where the 
majority of molecules are dissociated by F UV photons (e.g. in 
the nearby Orion Bar, iTielens et al.lll994l) . Hence, PAHs are 
a tracer of the atomic gas (for instance in the NG C 7023 neb- 
ula th ey follow the far infrared emission of C+, iJobUn et alJ 
I2OIOI) . and in the particular case of proplyds, PAHs will trace 
the morphology of the photo-evaporating flow. This is well 
illustrated in Fig. 3 showing the ionized (Ha), neutral (PAHs) 
and molecular (H2) gas components of the proplyd HSTIO. 



4.2. PAH abundance 

Proplyd models consider disk photoevaporation to be due 
mainly to efficient heating of the gas by energetic photo- 
electrons provided by small grains and PAH molecules. How- 
ever, these models assume PAHs in proplyds are as abundant 
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Table 1 

Parameters of NGC 7023 and HSTIO used to derive the PAH abundance in 
equation 1 



Parameter 




Reference 


NGC 7023 


^8.6 
A 1.25 


1306±5MJysr-' 
2723 ±6MJysr-' 


this paper 
this paper 


Go 


2.6 X 103 


Joblin etal. 2010 


Nh 


1 X 10^2 j.,„-2 


Joblin et al. 2010 


fPAH 

Jc 


7 X 10"^ 


Berne & Tielens 2012 


HST 10 


h.6 

hi .25 


720 ± 111 MJysr"' 
1466 ± 111 MJysr-' 


this paper 
this paper 




0.5 


this paper 


Go 


2.4 X 10^ 


Storzer & Hollenbach 1999 


Nh 


5.5 X lO"' cm-^ 


Storzer & Hollenbach 1999 


fPAH 

Jc 


8 X ur* 


this paper 



as in the interstellar medium (ISM), an assumption which has 
not yet been verified with observations. PAHs have been de- 
tected in, at most, 15% of the observed disks around isolated 
TTauri stars and they are u nder-abundant by a factor of 25 
when compared to the ISM (iGeers et al.ll2007l : lOliveira et al.l 
I20T0I) . The VISIR images of HSTIO presented in this paper 
can be used to estimate the abundance of PAHs in a proplyd 
PDR for the first time. 

The mid-IR PAH emission at a given wavelength is pro- 
portional to the num ber of carbon atom s locked in PAHs in 
the line of sight (e.g. iJobUn et"ani2010h . and on the intensity 
of the UV radiation field, GfQ Hence, we can write 



^ fPAH [^] »7 ^ 



(1) 



where /^'^'^ is the fraction of elemental carbon locked in 
PAHs, Nh is the column density of hydrogen atoms in the 
line of sight, [C]/[//] is the abundance of carbon relative 
to hydrogen atoms (1.6 x 10""^), and is the PAH emis- 
sivity at the given wavelength. The latter parameter, can 
be derived for sources where all the other parameters given 
in Eq. [1] can be determined independently as in the case 
of the reflection nebula NGC 7023. The values adopted for 
f^'^^, Nh and Go were taken from the literature (Table 1) 
while the intensities /g g and /11.25 for N GC 7023 were mea - 
sured directly in the Spitzer IRS spectrum (iPilleri et al.l2012l) . 
This yields for the PAH emissivities eg g = 4.5 x 10"'^ and 
fii.25 = 9.4 X 10""^ MJy sr^'cm^^Gy' which, when inserted 
in Eq. [T] combined with the parameters in Table 1, give the 
PAH abundance in the PDR of HSTIO, f^"^", if a good esti- 
mate of the column density of atomic gas in the line of sight 
Nh is provided. This parameter can be estimated from the 
electron density at the ion ization front which can be mea- 
sured with a fair accuracy. iBallv et al. (1998) obtain a value 
of «(, = 8 X 10^ cm"^ from the Y\a surface brightness, while 
IStorzer & Hollenbachl(ll999h find = 1.1 x 10^ cm ^ when 

' expressed in terms of the Habing field which corresponds to an integrated 
intensity between 9 1 .2 and 240nm of 1.6x10"' ergscm"^ s"' jHabind 19681) . 



correcting for the extinction to the Orion Nebula. Assum- 
ing pressure equilibrium at the ionization front, the density 
of H atoms in the neutral flow «//, must be of the order of 
10 X «(. or «H ~ 10* cm"^. For a head width of ~ 6 x 10'^ cm, 
measured in the PAH2 image, and assuming this value for 
the neutral flow lenght along the line of sight (symmetry in 
HSTIO), we obtain Nh ~ 6 x 10^' cm~^. From models and 
using the same parameters. fStorzer & Hollenbachl(ll999b find 
Nh — 5.5 X 10 cm"^. Adopting their value for Nh, and us- 
ing the brightnesses /g 6 and /11.25 (Table 1) measured in the 
VISIR PAHl and PAH2 images, we find similar PAH abun- 
dances in HSTIO of /^-^^ = 8.1 X 10""* and /^'^^ = 8.0x10"'*. 
These values are nevertheless extremely low: they correspond 
to an abundance of PA Hs 90 times lower than in NGC 7023 
jBerne & Tielensl2012l) . or about 50 times less than the value s 
found in the Orion Bar or in the diff^use ISM (lTielensl |2005|). 
And, since the brightnesses have been measured for the 
brightest pixels at the proplyd head (Sect. |3), our estimation 
gives the maximum abundance of PAHs in the evaporating 
flow toward 0' Ori C. 

5. DISCUSSION 

As other studie s of TTauri stars (iGeers et al.l 120071 
lOliveira et all 1201 Oh we find PAHs to be under-abundant in 
the PDR of the proplyd HSTIO, by a factor of 50 or more rel- 
ative to the diffuse ISM. The origin of this under-abundance 
cannot be readily explained, but some proposed hypotheses 
include clustering of PAHs followed by sedimentation inside 
the disk, destruction by FUV photons in the PDR, or de- 
struction by X-rays emitted by the low-mass central staiFl 
Nevertheless, more important than the causes for PAH under- 
abundance are the implications this result has on our under- 
standing of the physical processes shaping morphology and 
driving mass-loss in proplyds. Photo-evaporation of proplyds 
has been explained so far by the photo-electric heating of the 
gas which is known to have a much reduced photo-electric 
efficiency for grains larger than lOOA (lTielensll2005l) . Given 
the high column density Nh = 5.5 x 10^' cm"^ and low ex- 
tinction Av = 0.1 - 0.2 in the PDR of HSTIO, we expect 
these large grains to be depleted in the neutral flow. In fact, 
the photoelectric heating has it highest efficiency for PAHs 
(molecules of a few A) and the small end of the very small 
rains (up to a few tens of A). But according to IPilleri et al.l 
2OI2I) (Fig. 6 in their paper), for high radiation fields as those 
found in proplyds (> lO'^Go), the very small grains are de- 
stroyed and evaporated into free-flying PAHs. Therefore, we 
do expect PAHs to be the main agents of the photo-electric 
heating in proplyds. 

Assuming the disk surface has the same PAH abundance 
as the proplyd PDR (they are lifted from the disk surface by 
the evaporative wind), the low fraction of PAHs found in this 
letter will have a profound impact on the disk surface gas tem- 
perature. And hence, the current hypothesis of an evaporating 
disk creating the proplyd PDR and morphology may have to 
be revised. Two possibilities arise: 1) other gas heating mech- 
anisms are relevant for disk evaporation, such as collision al 
de-excitation of UV pumped H2 and H2 photo-dissociation 
followed by reformation on grain surfaces; 2) the PAH emis- 
sion from the proplyd cocoon is associated to remnant atomic 
gas from the protostellar envelope or the surrounding neb- 
ula. By combining detailed modeling with upcoming Her- 
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The majority of stars in the Orion Nebula cluster are in the mass range 
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of 0.1 - 1.5 Mo ([Hillenbrand & Hartmann 199: 
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schel (Berne et al., in prep.) and VLT data (Vicente et al., in 
prep.), the heating-cooling mechanisms at the disk surface can 
be assessed allowing to test each one of the possible scenarios 
creating the puzzling proplyd morphology. 
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